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Abstract 

Puzzling anomalies previously observed in multipulse NMR experiments in natural abundance ^^Si [A.E. Dementyev, D. Li, 
K. MacLean, and S.E. Barrett, Phys. Rev. B 68, 153302 (2003)] such as long-lived spin echoes and even-odd asymmetries, are 
also found in polycrystalline Ceo. Further experiments controlling the phases and tilting angles of the pulse trains, as well as 
analytical and numerical calculations allowed us to explain the origin of these anomalies. We prove that the observation of 
long magnetization tails requires two conditions: i) an rf field inhomogeneity able to produce different tilting angles in different 
sites of the sample and ii) the absence of spin diffusion (non-effective fiip-flop interactions). The last requirement is easily 
satisfied in diluted dipolar solids, where the frequency differences between sites, caused by disorder or other sources, are usually 
at least one order of magnitude larger than the dipolar couplings. Both conditions lead to the generation of stimulated echoes 
in Carr-Purcell (CP) and Carr-Purcell-Meiboom-Gill (CPMG) pulse trains. We show, both experimentally and theoretically, 
that the stimulated echoes interfere constructively or destructively with the normal (Hahn) echoes depending on the alternation 
or not of the n pulse phases in the CP and the CPMG sequences. Constructive interferences occur for the CP and CPMG 
sequences with and without phase alternation respectively, which are the cases where long magnetization tails are observed. 
Sequences with two, three and four vr pulses after the ■k/2 pulse allow us to disentangle the contributions of the different echoes 
and show how the stimulated echoes originate the even-odd asymmetry observed in both ^^Si and Ceo polycrystalline samples. 



During the last years, the control of quantum coher- 
ence has become of great importance for quantum in- 
formation processing and for the development of new 
nanotechnologies. In the nanoscale, many quantum phe- 
nomena become evident and can be exploited. However, 
interaction with various degrees of freedom of the envi- 
ronment degrades the quantum coherence. As a conse- 
quence, it is necessary to understand in detail the mech- 
anisms that produce decoherence and their timescales. 
Many suggestions to implement quantum computations 
involve nuclear or electron spins In particular, as sili- 
con has been an essential component of one of the main 
proposals,*' a group of researchers carried out a series 
of nuclear magnetic resonance (NMR) measurements to 
study the spin-spin "decoherence" time {T2) in ^^Si 
(4.67% natural abundance, spin i) doped with different 
concentrations of donors or acceptorssi Their results were 
completely unexpected and defied explanation to date. 

They performed T2 measurements with different pulse 
sequences commonly used with this purpose. They found 
that it is possible to detect ^^Si NMR signals up to much 
longer times (more than two orders of magnitude) by ap- 
plying the Carr-Purcell-Meiboom-Gill sequence as com- 
pared with the characteristic decay time obtained from 
the Hahn echo sequence ( « 5.6 ms). As these sequences 
usually give similar results, the difference cannot be dis- 
regarded and requires an explanation. 

Several tests were performed to determine the origin 
of this anomalous behavior. First of all, it was evident 
that there was no dependence on the amount of donors 
or acceptors in the ^^Si sample, even when this could 
change the linewidth and the spin-lattice relaxation time 
in one order of magnitude. Moreover, they found the 
same behavior in a pure single crystal of n.a. ^^Si?^ 

One hint was that the anomaly seemed to be present 
in systems with diluted spins. Thus, we implemented 



a series of related experiments and calculations in a 
Cgo sample. As each molecule has sixty carbon atoms 
and the magnetic isotope ^'^C, has a natural abundance 
1, 1%, in a typical sample there are 34.5% of molecules 
containing one ^^C, 11.4% containing exactly two, and 
2.4% containing three. 

As a first step, we search for an approximation to T2 
using the Hahn echo sequence^' [HE : (j) x ~ ''' ~ i'^)Y ~ 
T — echo], where r is a variable time. We obtained ~ 
15 ms. Then, we used different train pulse sequences to 
measure T2, derived from the Carr-Purcell sequence^ 
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and the Carr-PurccU-Meiboom-Gill sequenceji 
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Here, the numbers 1 and 2 have been added to the 
names of the sequences to denote one or two different 
phases for the tt pulses. 

We found that by applying the CPMGl sequence or 
the CP2 sequence, a magnetization tail appears, i.e. the 
signal remains for times much longer than T^^. More- 
over, with both of them we found that by choosing 2r 
-> j^HE ^]^g magnetization shows an even- odd behavior as 
observed in ^^Si^ (that is when odd echoes are already 
in the noise, even echoes are still visible). 

Another unexpected observation was the 
presence of noticeable stimulated echoes 
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when applying 
SE 



pulse sequences of the form, 

variable r and ti. In principle, these sequences involving 
TT pulses should not produce stimulated echoes. As we 
will show, the constructive or destructive interferences 
between the stimulated and normal echoes enable a phe- 
nomenological explanation for the long magnetization 
tails. 

Figure ^ shows the normalized magnetizations ac- 
quired by applying the sequences HE] CPl; CP2; 
CPMGl and CPMG2 to polycrystaUine Cqq. In the 
last four sequences t = 1ms and signal acquisition is per- 
formed at the top of the echo. In all the sequences the 
durations of the tt and it/ 2 pulses were carefully set from 
nutation experiments to 7.0 iis and 3.4 fis, respectively. 
No effort was done to reduce the rf field inhomogeneity 
which, for the completely full sample holder, was in the 
range 15-20%. It is clear from Fig.J^I that the magneti- 
zation, measured with the CP2 and CPMGl sequences, 
lasts for very long times («ls). 

On the other hand, by applying the CPl or the 
GPMG2 sequences the decay time is not longer than 
T^^, but it is remarkable that at short times the mag- 
netization makes a zigzag that, as will become clear, it is 
also a consequence of the presence of stimulated echoes. 
Besides, it should be noticed that the magnetization goes 
through negative values before reaching its final asymp- 
totic zero value. 
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FIG. 1: Comparison of the "C magnetizations acquired with 
different T2 sequences {t — 1 ms). The insert shows the 
decay observed with the Hahn echo, the decay simulated by 
considering only the Ising terms of the dipolar interaction and 
the short time regime of the CPMG2 sequence. 



The first question was if the T^^ measured for the 
Geo sample could be directly assigned to the dipolar in- 
teraction within the ^'^C nuclei. Our measurements of 



the spin-lattice relaxation time yielded Ti ~ 35 s reflect- 
ing that there are no efficient mechanisms for this kind 
of relaxation that should be taken into account in our 
analysis of the " decoherence" time. Thus, the main in- 
teractions in this system led us to write the Hamiltonian 
in the rotating frame as 
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where SuJi represents the chemical shift for spin i (relative 
to on-resonance spins); — )(1 — 3cos^6'y); 



/4 with 6*, 



vector r.j ■ and the static magnetic field 



the angle between the internuclear 
The gyromag- 
netic ratio of ^^C is 7 = 6.7283 x lO^rad T^^s^^ and the 
center to center distance between neighboring molecules 
is rij = 10 A.^ It should be noticed that because of the 
rapid isotropic rotations of the buckyballs in the solid, 
the spins can be considered to lie at the center of the 
molecule»ifi 

The Lorentzian full width at half maximum, 
FWHM ~ 167 Hz, corresponding to a free induction 
decay time Tj* ~ 1.91 ms is almost an order of magni- 
tude shorter than — 15.0 ms, thus reflecting the 
sample inhomogeneity. Then, it is usual to make the un- 



like spins approximation 



^ \5L0i — 5ujj\, and trun- 



cate the hij terms, leaving only the Ising contribution. 
Under these conditions, one can use the product oper- 
ator formalism, to obtain an analytical expression for 
the magnetization. Starting from the initial equilibrium 
density operator in the conventional strong field and high 
temperature approximations: 



p(i = 0)(x^/f 



(4) 



the application of the Hahn echo sequence with variable 
inter-pulse time r — kr, leads to the expression 
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(/y(fc2T)) = (0) m^o^K-^^) 



(5) 



Our calculation for a polycrystaUine sample of Ceo in- 
volves a set of 600 pieces 4000 unit cells each) with 
different representative orientations with respect to the 
external magnetic field and a random distribution of 
molecules with one, two or three spins each consider- 
ing the proportions expressed above. The characteristic 
decay time of the analytic curve is in very good agree- 
ment with the experimental one obtained from the Hahn 
echo sequence (see Ising model in Fig^ , manifesting that 
T^^ can be directly associated with the dipolar interac- 
tion. 

In order to look for the even-odd behavior reported for 
the Si sample^ in Ceo, we took long inter-pulse separa- 
tions (2t > T^^). We observed the anomalous behavior 
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FIG. 2: Even-odd effect observed with the CPMGl and CP2 
sequences in a Ceo polycrystalline sample. The interpulse 
separation is r = 8 ms. 
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FIG. 3: Stimulated echoes observed with the CPMGl^'^ and 
CPMG2^^ sequences in a Ceo polycrystalline sample. The 
parameters are r = 1 ms and ti — 15 ms, giving rise to a stim- 
ulated and a normal echo at 1 ms and 14 ms after the last 
pulse, respectively. The circles and squares represent splines 
(average of 20 original datapoints per symbol) performed to 
guide the eye. Notice the inversion of the phase in the stim- 
ulated echo. 



by applying the CPMGl and the CP2 sequences, as 
displayed in Fig|2] for r = 8 ms. While the first echo, 
occurring at t = 16 ms is higher than the second one, 
occurring at t = 32 ms, the third is lower than the forth, 
and so on. Moreover, the same sequence with t — 25 ms 
produces a first echo which is smaller than the second 
one. 

So far, all the results reported in silicon were also man- 
ifested in the Cgo sample, but there was still one more 
thing to check, the presence of stimulated echoes. 

We applied a sequence of three pulses^ 
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acq, 
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taking Lpi — X in all the experiments and varying 
and (^93 in correspondence with the sequences we used 
in the first experiments. For these stimulated echo se- 



quences, we denote \CPMG1 



SE 



^2 



[CP1^^:^2=^;^3 = ^], 
[CPMG2'^^ ■.if>2 = y\Vz = -y\ , and 
\(JP2SE : = X; = ~X\ . At a time ti - t 
after the third pulse, the normal echo (i.e. the refocusing 
of the Hahn echo at time t after the second pulse) 
appeared, and for all the sequences a stimulated echo 
peaked at time r after the third pulse. 

A very remarkable fact is that by applying the 
CPMGl^^ or the CP2^^ sequence, the stimulated and 
normal echoes have the same phases, but when applying 
the CPMG2^^ or the CPl^^ sequences they have oppo- 
site phases, as shown in FigOlfor the CPMG variations. 

Using the superoperator formalismjiSii^ 
for an arbitrary three-pulse sequence 



- T - {^2)^^ -h- (/?3)^3 - acgj where A 
and ifi are arbitrary tilting angles and phases in the 
XY plane respectively, we analytically solved the 
dynamical equations in the delta-pulse approximation. 
By taking into account that the FWHM ~ 167 
Hz and {■kT^^)~'^ ~ 21.2 Hz, we can neglect the 
spin-spin interactions and work with the Hamiltonian 

and the initial condition given in Eq.Q). Under these 
approximations, the following response function is ob- 
tained 
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where G{t) — &xp{itSujj) represents the normal- 

ized FID, the time t is the total time starting at the very 
first pulse, and the signal amplitudes fj = fjiPi, (32, Ps) 
are determined from the exact calculation: 
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The stimulated echo signal is the one that peaks at 
total time t = ti + 2t. Its amplitude is given by the ex- 
pression /e, evidencing that for sequences like the ones 
applied here (i.e. (32 = P3 = n), it cancels out. Conse- 
quently, the stimulated echo should not appear. 

The normal echo, second term in Eq.(|7I), peaks at a 
total time t ^ 2ti. Its phase in the CPMGl {ipi ^ 0, 
ip2 = 7r/2, Lpz = t:/2) or in the CP2 (i^i = 0, (^92 = 0, 
(y93 = tt) sequence is e*(2'^3-2v2+¥'i) = \ coinciding with 
that of the stimulated echo e*(^+¥'3+'^2-vi) ^ 1. However, 
if we apply the CPMG2 {ipi = 0, (^2 = 7r/2, tp^ = -it/2) 
or the CPl {ipi = 0, (p2 = 0, ip^ = 0) sequences they are 
in opposite phases, because the phase of the stimulated 
echo is inverted. Here, we have used ipi = for the X 
axis and it should be noticed that under the experimental 
conditions (/32 = P3) f2 and /g are always positive. Then, 
the remarkable fact is that these analytical calculations 
agree with the experimental result for any pair of tilting 
angles (32 = ^ tt. This lead us to conclude that not 
all the spins are affected by it pulses. We performed the 
three-pulse experiment for ipi = 1^2 = it /2, ip^ — tt /2, 
(3i = tt/2 and /?2 = f3z varying in the range [0, 57r/2] 
with T = 14 ms and ti = 15 ms. While the experiments 
reproduce the functional form of /g, there are no zeroes 
in the magnitude of the stimulated echoes for any value 
of P2- It is noticeable that the minimum magnitude of 
the SE, obtained for /32 = /Ja « tt is approximately 30% 
of the maximum value that occurs for P2 — Ps = tt/2 
(data not shown). 

It is evident that if we take ti = 2r in the three-pulse 
sequences © the normal and the stimulated echoes will 
appear at the same time, and as can be seeing in fig ||2J), 
they will interfere constructively in the CPMGl and the 
GP2 sequences or destructively in the CPMG2 and the 
GPl ones. As a consequence of these constructive inter- 
ferences and the well known fact that stimulated echoes 
have much longer decay timesjSiiS the magnetization tails 
appear. 

In order to analyze the different contributions to the 
even-odd asymmetry in the echoes, displayed in Fig[21 we 
performed experiments with two, three and four tt pulses 
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FIG. 4: Echoes obtained after three and four tt pulses with 
the (a) CP2^-^ and (b) CPMCl-^-^-like sequences. The time 
between the tt/2 and the tt pulse is r = 3 ms, while the 
intervals between tt pulses are ti = 10 ms. For three tt pulses 
the stimulated echo appears at ti — r = 7 ms after the last 
pulse, while the normal echo appears at r. The normal and 
stimulated echo times after the last tt pulse are exchanged 
in the case of four tt pulses. Notice that after four tt pulses, 
the normal echoes have almost dissapeared. The circles and 
squares represent splines (20 original datapoints averaged). 



after the 7r/2 pulse with interpulse separations preclud- 
ing the overlap of the normal and stimulated echoes. In 
particular, we set the separation between the 7r/2 and 
the TT pulses r = 3 ms, while the separations between tt 
pulses where ti = 10 ms. The results for three and four 
TT pulses for the case of a CP2^^ and CPMGl^-^-like 
sequences are shown in Fig. 0| (a) and (b) respectively. 

As can be seen, the stimulated echo occurring after 
four TT pulses, i.e. at total time t — 3ti -I- 2r or r af- 
ter the last pulse, is larger than the one after three tt 
pulses t = 3ti(or ti — t after the last pulse) in both 
cases. We use the same formalism and Hamiltonian ap- 
plied for the three-pulse sequences, setting a single tilting 
angle (3 for all the pulses after the 7r/2 one, with (pi — 
and Pi — or tt alternatively for z 7^ 1 {GP2^^) or 
p>i^Q and p>, = 7r/2 for i ^ 1( CPMG\^^). Thus, we 
obtained for the stimulated echo amplitudes the func- 
tions: ^47r = l/2sin2/3(5/2cos2/3 - cos/3 -h 1/2) and 
Astt = =Fl/2sin^ /?(1 — cos/3). The minus sign in As^r 
corresponds to the CP2^^ sequence. Again, it is evident 
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that for /3 = TT, there should be no stimulated echo, i.e. 
AiTT = ^37r = 0. However, it can be seen from these for- 
mulae that for angles in the interval f3 = 7r(l ± 0.35), one 
gets AiT^ > ^3^, which is in qualitative agreement with 
the experimental observation. Numerical calculations for 
longer sequences (i.e. a larger number of tt pulses) sup- 
port this alternation of amplitudes. Consequently, hav- 
ing proved that the stimulated echoes interfere construc- 
tively with the normal echoes in the sequences CPMGl 
and CP2, we have found the explanation for the even- 
odd asymmetry observed in Fig|21 

In order to further verify the fact that not all the sites 
in the sample are experiencing the same externally fixed 
tilting angle, we applied the same sequences and (0) 
but applying our best 2tt or An pulses (zeroes in the nuta- 
tion experiment), instead of tt pulses. In all of them, we 
kept the first {■§) ^ pulse. Under these conditions, if the 
pulses were perfect, one would expect to acquire just an 
FID signal. This is exactly what we got by applying the 
CPMG2(^2-k) or the CPl(2,r) sequences or their equivalent 
with 47r pulses. Instead, by applying the CPMG1(^2-k) 
or the CP2(27r) sequences the magnetization tails reap- 
peared as also occurred with the An pulses. 

The conclusion that not all the spins experience the 
same tilting pulse should not be surprising, as it is well 
known that the rf field inhomogeneity may be around 5- 
10% in ordinary experimental conditions. However, the 
remaining question is why these long tails do not appear 
when working with magnetically concentrated samples,^"* 
as usually occurs in solid-state -'^H-NMR measurements. 
While the rf field inhomogeneities occur in all cases, in di- 
luted systems, the dynamical contribution of the dipolar 
interaction (i.e. the flip-fiop part) is not effective. Then, 
it can be conceived that in concentrated systems, there 
are different tilting angles in different sites, but during 
the interpulse time the defects or excesses in the tt pulses 
are exchanged among the sites because of the flip-fiop, 
thus, effectively averaging to the externally fixed pulse. 
By contrast, in diluted systems the imperfection over one 
site has always the same sign or direction, originating the 
stimulated echoes and consequently the long tails. This 
fact was verified through numerical simulations where 
imperfections up to 10% in the "tt pulses" were ran- 
domly applied to different spin sites. But a very impor- 
tant ingredient in this explanation is that the disorder in 
energies among the different sites must be greater than 
the dipolar couplings. This condition occurs quite nat- 
urally in diluted systems because the dipolar interaction 



is very small; thus, any interaction capable to produce 
differences in site energies will be enough to hinder the 
flip-flop dynamics. As a practical warning, one should 
consider the fact that magnetically diluted samples re- 
quire larger amounts making more serious the problem 
of rf field inhomogeneity. 

The long tails obtained in our experiments can not be 
interpreted as a decrease in the decoherence rate of the 
spin system. Because of the pulse imperfections, some 
magnetization is preserved in the direction of the external 
magnetic field, having a decay time in the order of the 
spin-lattice relaxation time>i^ 

Our findings may be significant for qubit coherence 
control in a nuclear spin bath, as occurs with electron 
spins in a lattice containing Si, GaAs, In or Al.^^-^^ In 
these situations, the electron spin coherence is limited 
by spectral diffusion arising from dipolar fiuctuations of 
lattice nuclear spins. 

Another important aspect to discuss is that, as it is 
well known, the dipolar dynamics is reversible ^^iiSiiS and 
there is no real loss of information in the collective spin 
system because of it. In this sense, a more relevant mea- 
surement of "decoherence time" in Cgo, would involve 
an homonuclcar decoupling as performed by Ladd et al^ 
in silicon. Under those conditions, more detailed experi- 
ments are required to test if there is still a contribution 
from stimulated echoes. We find no reason to discard it 
beforehand, because the necessary conditions: rf inho- 
mogeneity and no ffip-fiop dynamics will still be present. 

The absence of nuclear spin diffusion observed in 
our experiments and in those performed in n.a. ^^Si^ 
seems to be a particular manifestation of the localization 
phenomenon^SiSi However, as one is working with a 3-d 
many-body system, the complexity of a theoretical ap- 
proach is beyond the scope of this work. Nevertheless, 
NMR studies as a function of nuclear spin concentrations 
and using different time windows may help to search for 
a possible phase transition, opening up a new way to look 
at localization phenomena. 
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